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We analyze two recently proposed measures of non-Markovianity: one based on the concept of 
divisibility of the dynamical map and the other one based on distinguishability of quantum states. 
We provide a toy model to show that these two measures need not agree. In addition, we discuss 
possible generalizations and intricate relations between these measures. 
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I. INTRODUCTION 

The dynamics of open quantum systems attracts nowa- 
days increasing attention It is relevant not only 
for the better understanding of quantum theory but it 
is fundamental in various modern applications of quan- 
tum mechanics. Since the system-environment interac- 
tion causes dissipation, decay and decoherence it is clear 
that dynamic of open systems is fundamental in mod- 
ern quantum technologies, such as quantum communi- 
cation, cryptography and computation Q. The usual 
approach to the dynamics of an open quantum system 
consists in applying an appropriate Born-Markov approx- 
imation leading to the celebrated quantum Markov semi- 
group [fl 0] which neglects all memory effects. However, 
recent theoretical studies and technological progress call 
for more refine approach based on non-Markovian evolu- 
tion. 

Non-Markovian systems appear in many branches of 
physics, such as quantum optics [l|, 0] , solid state physics 
[8(, quantum chemistry [9j, and quantum information 
processing (Io| . Since non-Markovian dynamics modifies 
monotonic decay of quantum coherence it turns out that 
when applied to composite systems it may protect quan- 
tum entanglement for longer time than standard Marko- 
vian evolution In particular it may protect the sys- 
tem against the sudden death of entanglement [T3|. It 
is therefore not surprising that non-Markovian dynamics 
was intensively studied during last years [l3j . 

Surprisingly, it turns out the concept of 
(non)Markovianity is not uniquely defined. One 
approach is based on the idea of the composition law 
which is essentially equivalent to the idea of divisibility 
[lH ]. This approach was used recently by Rivas, Huelga 
and Plenio (RHP) |l6| to construct the corresponding 
measure of non-Markovianity, that is, RHP measure 
the deviation from divisibility. A different approach 
is advocated by Breuer, Laine and Piilo (BLP) in 
Ref. BLP define non-Markovian dynamics as a 

time evolution for the open system characterized by a 
temporary flow of information from the environment 



back into the system and manifests itself as an increase 
in the distinguishability of pairs of evolving quantum 
states. It is clear that RHP characterize a mathematical 
property of the dynamical map whereas the idea of BLP 
is based on physical features of the system-reservoir 
interaction, rather than the mathematical properties of 
the dynamical map of the open system. In a recent paper 
[ri| Haikka, Cresser and Maniscalco performed detailed 
analysis of these approaches studying the dynamics of 
the driven qubit in a structured environment. It is 
indicated [17| that the concepts of RHP and BLP need 
not agree, as was also conjectured in [l8| and checked 
in [l9(. BLP measure was recently analyzed for the 
dynamics of a qubit coupled to a spin environment via 
an energy-exchange mechanism (20| . 

In the present paper we perform further analysis of this 
problem. In particular we provide a simple toy model 
showing that Markovian evolution a la BLP may be in- 
divisible and hence non-Markovian according to RHP. 
Actually, one may feel the relation BLP vs. RHP quite 
analog to the relation between separable and PPT states 
in entanglement theory, where separable states define the 
proper subset of PPT states. States which are PPT 
but entangled are bound entangled. Using this anal- 
ogy one may call Markovian evolution a la BLP but in- 
divisible (non-Markovian according to RHP) - 'bound 
non-Markovian'. Finally we discuss possible general- 
izations of BLP using tensor product structures, and 
show that a slight modification of BLP reduces their 
concept of Markovianity to divisibility (the bound non- 
Markovianity can be washed out). The main motivations 
of this paper are to expose the relation between two dif- 
ferent concepts on a simple example, and to present the 
way to unify them, which we hope clarify the issue. 



II. DIVISIBILITY VS. BACKFLOW OF 
INFORMATION 

The measure for non-Markovianity proposed by Rivas, 
Huelga and Plenio (RHP) [l6| is based on a notion of 
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divisibility: a trace praserving completely positive map 
A(f, 0) is divisible if it can be written as a 



A(t + T,0) = A(f + r,f)A(f,0) 



(1) 



and A(f + r, f ) is completely positive for any t, r > 0. 
RHP define a map to be Markovian exactly when it is 
divisible. Note that 

A(t + r,t) = A{t + T,0)A-\t 1 0) , (2) 

satisfies composition law 

A(s,t) = A(s,u)A(u,t) , (3) 

for any s > u > t, which is usually attributed for Marko- 
vian evolution. It is shown fl6l ] that the quantity 



7(f) = lim 

e->-0+ 



|ll d ®A(f + e,f)P/||i-l 



(4) 



enjoys g(t) > if and only if the original map A(f, 0) is 
indivisible ( 1 1 • 1 1 1 is a trace-norm) . As usual Id denotes an 
identity map in Md(C), and denotes the maximally 
entangled state in C d ® <C d (we consider the evolution of 
a d- level system). Now, using the fact that any divisible 
(and differcntiable) completely positive map satisfies a 
local in time master equation 



-A(f,0) = L(f)A(f,0) 



A(0,0) = ILj 



(5) 



where the local generator L{t) is a legitimate Markovian 
generator for any t > 0. The formula Q may be equiv- 
alently rewritten in terms of L(t): 



g(t) 



lim 

E^0 + 



l d ®(l d -eL(f))]P+|| 1 



(6) 



A second criterion of non-Markovianity was proposed 
by Breuer, Laine and Piilo (BLP) in Ref. The 
BLP criterion identifies non-Markovian dynamics with 
certain physical features of the system-reservoir interac- 
tion. They define non-Markovian dynamics as a time 
evolution for the open system characterized by a tempo- 
rary flow of information from the environment back into 
the system. This backfiow of information may manifest 
itself as an increase in the distinguishability of pairs of 
evolving quantum states. Hence, according to BLP the 
dynamical map A(f , 0) is non-Markovian if there exists a 
pair of initial states p\ and p 2 such that for some time 
t > the distinguishability of p\ and P2 increases, that 
is, 



a(pi,p 2 ;t) = —D{p 1 (t),p 2 (t)] 



(7) 



where D(pi,p 2 ) = — P2II1 is the distinguishability 

of pi and p 2 , and p k (t) = A(t,Q)p k . 



Using these two criteria one easily defines the corre- 
sponding non-Markovianity measures of the dynamical 
map A: 



A/rhp (A) = 



1+1 ' 



where X = g(t)dt. Similarly, one has 



A/blp(A) = sup 



Pl,P2 J a>0 



a(p 1 ,p 2 ;t)dt 



(8) 



(9) 



It is clear [15| that A/rhp (A) = implies A/blp(A) = 0. 
Hence, all divisible maps are Markovian according to 
BLP. The converse is in general not true. This prob- 
lem was carefully analyzed in the paper [ITjj for phc- 
nomenological integro-differential master equations and 
in the recent paper [i~7j . where the authors have stud- 
ied the non-Markovian character of a driven qubit in a 
structured reservoir in different dynamical regime using 
A/rhp (A) and Ablp(A). 



III. TOY MODEL - CLASSICAL STOCHASTIC 
DYNAMICS 

In should be clear that the above discussion applies for 
a classical stochastic systems as well. Consider stochas- 
tic dynamics of d-state classical system described by 
the probability vector p = (j>i, . . . ,p^). Its evolution 
p(0) — > p(f) = A(f) ■ p(0) is defined by the family 
of stochastic matrices A(f), that is, A.y(£) > for all 
i,j=l,...,d and 



(10) 



1=1 



for each j = 1, . . . , d and t > 0. These conditions guar- 
antee that p(f) is a legitimate probability vector for all 
t > 0. The corresponding local-in-time master equation 



Pi(f) =Y j L ij {t) Pj {t) 



(11) 



is defined in terms of the local generator L(t) which is 
defined in by the classical dynamical map A(t) as follows 



L(t) = A(f) • A-\t) 



(12) 



In particular if L does not depend on time A(f) defines a 
one-parameter semigroup of stochastic matrices A(t) :~ 
A(t, 0) = e tL . Let us recall that A is stochastic if if and 
only if it satisfies well known Kolmogorov conditions: 



Lij > 



i ^ j , and ^ L v} = 



(13) 
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for each j = 1, . . . , d. Now, a stochastic evolution A(t, 0) 
is divisible if it can be written as a composition of two 
stochastic maps A(t + r, t) and the original A(t, 0) 



A(t + T,0) = A(t + T,f)-A(t,0) , 



(14) 



for any r > 0. Clearly, A(i, 0) is divisible if and only if the 
corresponding local generator L{t) satisfies Kolmogorov 
conditions for any t > 0. Now, we would like to compare 
two non-Markovianity measures due to RHP and BLP. 

It is clear that classical evolution may be rewritten 
in the quantum framework as follows: any probabil- 
ity vector p gives rise to the diagonal density matrix 
p = 'Yl,kPk\k){k\ 1 and hence the stochastic map - clas- 
sical channel - A gives rise to the following Kraus repre- 
sentation 



i,j=i 



(15) 



Now, to apply ((4]) one needs the classical analog of 
(l<g)A)P + . Let us define (H cl <g> A)P C |, where the clas- 
sical identity map is defined by 



W = E \ i )( i \p\ i )( i \ > 



(16) 



and the classical analog of the maximally entangled states 
reads as follows 



One finds therefore 



1 



(]l cl ®A)P+ = - AjimiSliM 
Now, the 'classical' quantity 



(17) 



(18) 



hm IUci®A(t + M)iSlli-l 

e^0+ e 



is strictly positive if and only if the original stochastic 
map A(i,0) is indivisible and hence g(t) identifies non- 
Markovianity of the classical stochastic process. Equiva- 
lently, it may be rewritten using local generator L(t) 



g(t) = lim 



|[l cl ®(]l cl - eL(t))]P c 



ilH 



1 



(20) 



Consider now the following toy model of stochastic dy- 
namics of 2-state system 



A(t,0) 



l-Xo(t) Xi(t) 
X (t) l-Xi(t) 



(21) 



where xo(t),xi(t) G [0, 1] for all t > 0. One finds for the 
local generator 



L(t) = 



-a (t) ax(t) 
a (t) -ai(t) 



(22) 



where 



(t(. 



±p(l - Xi) + XjXp 
1 — Xq — X\ 

AqXi + ±i(l - Xp) 
1 — Xq — X\ 



(23) 
(24) 



Now, A(t, 0) is divisible if and only if ao(t), ai(t) > for 
alH > 0. On the other hand one easily finds 



a(t) = -2[oo(t) + oi(t)]|A | , 



(25) 



where Aj. = (pi — P2)k- Hence, contrary to g(t), the 
quantity a(t) controls only the sum of eto and a\. In 
principle one may have a(t) < even if ao(t) < or 
ai (i) < 0. Indeed, let us define 



f ( T )dr , Xl (t) = / h{r)dT , (26) 



x (t) -- 

such that < J Q /fc(r)rfr < 1 , for all t > 0. Hence, our 
toy model is fully controlled by functions /o and fi . Let 

fo(t) = Ksini , t>0, (27) 

and /i(f) = for t G [0,7r] together with 

= -Ksini , i > 7T , (28) 

where < k < 1/2. One finds 

K sini 



OoW + ai(0 



1 — K + ft cos £ 



*G[0,7r], (29) 



and ao(t) + ai(t) = for t > n. Note, that for t > ir one 
has 



a o(t) = —a-i(t) = ft sin t , 



(30) 



which proves that A(t, 0) is not divisible. However, this 
map is 'classically' Markovian according to BLP [15( due 
to o (i) + ai(t) > for all t > 0. 



IV. QUBIT DYNAMICS 

Consider now the following dynamics of a qubit 

poo(t) = pooX (t) +p n [1 - xi(t)] , 

pu(t) = poo [1 - xo(t)] +pnxi(t) , (31) 

poi(t) = Poil(t) , 

where xo(t),xi(t) G [0, 1], and 

\j{t)\ 2 <xo{t) Xl (t) . (32) 

The above conditions for Xk(t) and j(t) guarantee that 
the dynamics is completely positive. One easily finds for 
the corresponding local generator 

D, 1 y 

Lp = -i — [<T z ,p] + '^2a k L k p+-L z p, (33) 

k=0 
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where 



L p = a+pa- - -{cr_cr + ,p} 
1 r 

Lip = (T-P<J+ - -{(T + a-,p\ 
L z p = <J z pa z - p . 



(34) 



The time-dependent coefficients do and ai are defined in 
([23]) and (|24|) . whereas r(i) and f2(i) read as follows 



r(t) 



Re 



7(0 



n(t) = im^f 

7W 



(35) 
(36) 



The corresponding dynamical map A(t, 0) is divisible if 
and only if a k (t),T(t) > for all t > 0. On the other 
hand one finds the following formula 



a{t) 



2A{t)^ + [A{t)+AT{t)]\A m {< 



(37) 



VAg + |A i| 2 

where A(t) = ao(t) + ai (£) , and a 2 x 2 matrix A reads 

A = Pl (0) - P2 (0) . (38) 

It is therefore clear that using the same arguments as in 
the classical toy model we may have a(t) < but the 
dynamical map is not divisible. 

V. GENERALIZED CONDITIONS 

We have studied so far the non-Markovian character 
of simple classical and quantum (toy) models and com- 
pare two non-Markovianity measures due to RHP and 
BLP. We performed explicit construction showing that 
in gene ral these two measures do not agree supporting 

|i5Hig|. 

Let us observe that both approaches may be analyzed 
from a slightly more general perspective. Consider a fam- 
ily A(i, 0) of trace preserving positive maps (not necessar- 
ily completely positive). Actually, the approach of BLP 
needs only positivity of A(i,0). It is clear that A(t,0) is 
contractive (one has | |A(i, 0)p\ i = |H|i). Now, A(t, 0) 
is divisible if 

A(t + r,0) = A(t + r,f) ■ A(f,0) , (39) 

and A(t + r, t) is a contraction for all r > 0. It implies 

||A(t + r,i)|| = l , (40) 

where the norm of the map is defined via 

1 1 A|| = sup ||Ap||i . (41) 
IWIi=i 



Note, that A(t, 0) is divisible if and only if 
d 



dt 



A(t,0) = L(t)A(t,0) , 



(42) 



where L(t) is a local generator satisfying the following 
dissipativity condition [2l| (see also Q): for any set 
of mutually orthogonal rank-1 projectors Pj such that 
J2 i Pi = I, one has 



Tt[P<L(i)P 3 -] > , ijLj, 
^Tr[pi(i)P,] = 0. 



(43) 



The second condition guaranties that A(t, 0) is trace pre- 
serving. Hence, formula is equivalent to 



lim 

E^0 + 



\t-eL(t)\\-l 



= 



(44) 



We call A(t, 0) completely divisible iff A(t + t, t) is com- 
pletely contractive. Recall that $ : B(T-L) — > B(H) is 
completely contractive [22| iff life <£> $ is contractive for 
all k = 2,3,.... If dim'H = d < oo, then, as in the 
case of complete positivity, it is enough to check that 
Id <g) $ is contractive. Note, that due to the fact that 
A(t + t, t) is trace preserving complete contractivity is 
equivalent to complete positivity. Hence, complete di- 
visibility is the standard divisibility considered by RHP 
[l6| . Finally, recall that $ is completely bounded [22J iff 
% <g> $ is bounded all k = 2, 3, — One defines 



and 



1*1 1 (fc) = l|lfc®$ll , 



| $ | | cb = supply 

k 



(45) 



(46) 



$ is completely contractive iff ||$|| c b < !• If dim'H = d < 
oo, then necessarily $ is completely bounded. Moreover 



1*1 



cb 



1*1 



(d) 



Hence, A(t, 0) is completely divisible iff 
||A(t + T,t)|| cb = 1 , 



(47) 



(48) 



which is much more restrictive than (|40p . Interestingly, 
for completely positive * one has 



*llcb = ll(ld®*)Pj 



(49) 



which reproduces result of RHP. Actually, one may gen- 
eralize the whole analysis to (completely) contractive dy- 
namics on the Banach space generalizing old results for 
contractive semigroups [231 ] - 

Now, let us observe that the approach of BLP which is 
valid for positive maps A(i, 0) may be easily generalized 
if we assume that A(t, 0) is completely positive. One may 
introduce 



o-(pi,p 2 ;t) = — D[ Pl (t),p 2 (t)} 



where 



p fc (t) = [l d ®A(t,0)]p k 



1,2 . 



(50) 



(51) 
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It is clear that if a(pi,p2',t) < for all initial states 
Pk, then a(p±, p2',t) < for all initial states pk- Note, 
however, that the converse needs not be true. Again, 
one may show that there exists indivisible maps A(t, 0) 
for which a(pi,p2;t) < for all initial states pk- Note, 
that condition a(pi,p2',t) < may be reformulated as 



2F(A;t) 



dt 



|A(i)||i <0 



(52) 



where A(t) = [l d <g> A(t, 0)]A, and A t = A together with 
TrA = 0, which follows from A = (pi — f>2)/2 and pk are 
true states. It should be stressed that condition TrA = 
is very restrictive. Let us recall that if $ is trace preserv- 
ing then $ is positive iff 



\®a\U < \\a\U 



(53) 



for all a) = a. Note however that if one restricts to a 
enjoying Tra = 0, then ||$a||i < ||a||i does not imply 
positivity of $. If we relax TrA = 0, then the condi- 
tion a (A; t) < for all Hermitian A imply that A(t, 0) is 
(completely) divisible and hence definition of Markovian- 
ity due to BLP reduces to divisibility. More explicitly this 
fact can be formulated as: 

Theorem. For a bijective evolution A(t,0), a(A;t) < 
for all A^ = A, if and only if it is divisible. 

Proof. The 'if part is straightforward because the trace 
norm is monotonically decreasing given the completely 
contracting property of the A(t + r,t), 

||A(t + r)||i = || [Id ®A(t + T,t)]A(t)|[i < IIA^H: , 

for every t and r > 0. Conversely, if the evolution A(t, 0) 
is bijective the partitions A(i+r,t) = A(t + r, 0)A _1 (t, 0) 
are well-defined. Then ct(A; t) > for some t implies that 
for some small r 

||A(t + r)||i = ||[]l d ®A(t + x,i)]A(i)|| 1 > ||A(t)||x. 

So there exists some partition A(t + r, t) which is not 
completely contractive, and therefore the dynamics is not 
divisible. 

□ 



VI. OPERATIONAL MEANING OF 5(A;t) 

One may think that by relaxing the condition TrA = 
the interpretation of an increase of <r(A;i) in terms of 
backflow of information from environment to system is 
lost, however that is not the case. To explain this more 
carefully we have to recall here some results from quan- 
tum hypothesis testing, particularly the one-shot two- 
state discrimination problem [24| . 

Consider a quantum system whose state is represented 
by the density matrix p\ with probability p, and pi with 



probability (1 — p). We want to determine the density 
matrix that describes the true state of the quantum sys- 
tem by performing a measurement. If we consider some 
general POVM {Mj}, where j £ H is the set of possible 
outcomes, we may split this set in two cases. If the out- 
come of the measurement is inside of some subset Acfi, 
then we say that the state is p\ . Conversely if the result 
of the measurement belongs to the complementary set A° 
such that A U A c = fi, we say that the state is p2- Let us 
group the results of this measurement in another POVM 
given by the operator T = J2j£A Mj- 

Thus, when the true state is p\ (which happens with 
probability p) wc erroneously conclude that the state is 
P2 with probability 



(1-p) £ TrfaMj] = (l-p)Tr 

= (l-p)Tr[ Pl (I-T)]. 

On the other hand, when the true state is p2, we erro- 
neously conclude that the state is p\ with probability 




3&A 



f>2 



= pTr [p 2 T] 



Note that when p = or p = 1 wc immediately ob- 
tain zero probability to identify wrongly the true state. 
The problem in one-shot two-state discrimination is to 
examine the tradeoff between the two error probabilities 
pTr [p 2 T] and (1 - p)Tr [p x {l - T)]. Thus, consider the 
best choice of T that minimizes the total error probabil- 
ity 

^m^pTr [ P2 T] + (1 - p)Tr [ Pl (l - T)]} 

= ™^ I {( 1 ~P) + Tr \PPz T - C 1 - P)Pi T }} 
= (l-p)- max [Tr (AT)], 

0<T<I 

where A = (1 — p)pi — pp2 is a self-adjoint opera- 
tor (sometimes called Hclstrom matrix (25|) with trace 
TrA = 1 — 2p, which vanishes only for the unbiased case 
p = 1/2. By using analogous arguments that for p = 1/2 
(see [j, HH) the result of the optimization process turns 
out to be 

mm^pTr [p 2 T] + (1 - p)Tr [ Pl (l - T)}} = 1 ~ ^ . 

(54) 

Thus the trace norm of A = (1 — p)px — pp2 gives 
our capability to distinguish correctly between p\ and p2 
in the one-shot two-state discrimination problem. Since 
that is only a function of the information we gather by 
the prior probability p and by the measurement T, ||A||i 
is a measure of the information we have. 

Consider again a trace preserving map $ : B(H) — > 
B(H). If it increases the trace norm of A, ||$(A)||i > 
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|| A ||i, we can assume that $ carries information about 
the correct state of the system. Otherwise it cannot de- 
crease the probability to identify wrongly the true state!. 
For that reason, if we process data before making a mea- 
surement, which means to apply some 4> over the states, 
we cannot increase the trace norm of A; since we cannot 
gain more information about some data just by process- 
ing it!. This is the so-called data processing inequality 
[24j| . and implies that $ is contractive, i.e. positive. In 
addition, we can never dismiss that the states p\ and p 2 
are part of a larger system "SA" . In such a way that 
they are the result of tracing out the ancillary systems 
px,2 = ^ a[pia,2A\- By making the same analysis as be- 
fore, in order not to gain information we then require the 
map id (g> $ to be contractive, i.e. $ completely contrac- 
tive and thus completely positive. 

Going back to the problem of Markovianity, one may 
now realize that (t(A; t) is a measure of the information 
gained by the system for some initial Hclstrom matrix 
A. This interpretation puts together the two concepts 
of Markovianity based upon divisibility and backflows 
of information. Particularly if one ignores the potential 
existence of an ancilla, under the case of unbiased dis- 
crimination problem, ||A||i = ^||pi — P2H1, which is the 
definition of the trace distance between pi and p 2 , re- 
covering the approach of BLP. Nevertheless in principle 
there is no reason to presuppose the fulfillment of these 
particular assumptions. 

One can define a measure of non-Markovianity as ([9]) , 
with ct(A;<) in the place of o{p\ 1 p%\ t). That will be 
zero iff © is, however §5§ is easier to compute because 
it avoids the complicated optimization procedure in (|9]) 
now upon every Helstrom matrix A. 

We have already proposed examples of non-divisible 
dynamics with Ablp(A) = 0, but for completeness we 
have computed ct(A; t) for the phenomenological integro- 
differential master equation [26[: 



dp(t) 
dt 



where k(t) = 7 * and 



k(t')Lp(t - t')dt' , 



Lp = Jo{n + l)L p + "f nLip, 



(55) 



(56) 



with 7,7o,n > and L$ and L\ were given in (|34|) . De- 
spite this kind of equations is used as a simple model 
of non-Markovian dynamics, it was showed in [l9| that 
A/blp(A) = 0. However a(A;t) is not always decreasing 
for any Helstrom matrix. For example by taking p = 0.07 
and 




FIG. 1: Behavior of cr(A; t) for the integro-differential model 
(|55[) . The initial Helstrom matrix is A = 0.93/5i — 0.07/92, 
where pi and p2 are given in (|57l) and (|58[) . There is a time 
period between 100 and 150 where <t(A; t) increases. 



there is a period where (t(A; t) grows as we have repre- 
sented in Fig 1. This denotes the existence of a backflow 
of information which was expected from the non-divisible 
character of this dynamics. 



VII. CONCLUSIONS 



We have analyzed two concepts of Markovianity one 
based on the divisibility property of the dynamical map 
and the other based upon the distinguishability of quan- 
tum states. We have given very simple examples where 
these two criteria do not coincide. Furthermore we have 
proposed a way to make them equivalent, in the sense 
that Markovianity would be identified by divisibility, but 
keeping the interpretation in terms of flows of informa- 
tion. For that we resort to the results in the one-shot 
two-state discrimination problem; and point out that an 
increase of the trace norm of any Hcrmitian matrix dur- 
ing the dynamics can also be associated with a backflow 
of information from the environment to the system. 
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